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Abstract: The paper studies spin-orbit interaction (i.e. the effect the spin has on the particle's trajectory 
in a magnetic held) as a model of quantum computation. The two-level spin quantum system is examined 
using the stochastic mechanics formulation. The control of the entangled spin state is considered as a 
problem of control of the mean moment of a particles ensemble along a reference trajectory. It is shown 
that such a control can be succeeded by applying an open-loop control scheme. 

1 Introduction 

The spin states, i.e. the eigenstates associated with the measurement of a particle's magnetic moment along 
the z-axis are a fundamental element of highly entangled states. For example, in quantum computing all 
the unitary entangle operations on many spins can be implemented by compositions of those on the two 
spins (universality of quantum circuits). Therefore, the stabilization of the two-spin system is of great 
importance for designing a quantum entangler machine. There are two types of control for chaotic and 
quantum systems, i.e. open-loop control where the controls are predetermined at the start of the experi- 
ment and closed- loop (or feedback control) where the control can be chosen through-out the experiment [Tj . 
Previous work on quantum open- loop control includes flatness-based control on a single qubit gate [2] . In 
this paper the classical Schrodinger's equation is used and quaternion notation for the spin's Hamiltonian 
is applied. Then a flatness-based controller is derived for steering the particle's transition between the 
spin's eigenstates. On the other hand, quantum feedback control was developed as a quantum analogy 
to the classical theories of nonlinear (Stratonovich) filtering [3]. This approach is based on an analogous 
of the separation principle which holds in classical stochastic control [4]. Quantum feedback control is 
actually observer-based control. First quantum filtering is used to obtain an estimation of the stochastic 
quantum variable and then a feedback controller is designed based on the output of the quantum filter. 
The quantum filter describes a classical stochastic process [5] , [7] , [8] . 



The present work studies spin-orbit interaction as a model of quantum computation (i.e. the effect the spin 
has on the particle's trajectory in a magnetic field) and examines control of the two-level spin quantum 
system using the stochastic mechanics approach. In stochastic mechanics each particle follows a continuous 
path which is random but has a well-defined probability distribution [9] , [10] . Stochastic mechanics has been 
recognized as a self-consistent formulation of quantum mechanics in the framework of stochastic processes 
and is established as part of stochastic control theory [IT] , [12] , [13] , [14] . The particles' kinematic model 
is associated to the model of a quantum oscillator, taking the particle to be a 3-DOF variable (position 
in a cartesian coordinates frame) |15|.[16|. When the effect of the spin in the particle's motion is also 
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considered, and the particle's motion takes place under an external magnetic field (as described in the 
Stern-Gerlach experiment) then the deflection of the particles' trajectory should be also taken into account 
, Q~7] . The component of the spin along the field becomes a discrete random variable which is correlated 
with the average velocity [TT] . 

The concept adopted in this paper is that starting from Schrodinger's equation and passing through the 
Fokker-Planck equation and the Ornstein-Uhlenbeck diffusion one finally arrives at Langevin's stochastic 
differential equation. A form of Langevin's equation is also used to describe the variations between the 
spin-eigenstates (discrete energy levels of the spin system) [TT].[12j. This enables to consider the control 
of the entangled spin state as a problem of motion control of the mean of a particles ensemble along a 
reference trajectory. It is shown that using an appropriate open-loop control scheme (which is based on 
flatness-based control theory and which can be realized by a magnetic field) motion control of the particles 
mean along the desirable trajectory is possible. This in turn implies that the quantum state represented 
by the particles moment can be also controlled. 

2 The spin as a two-level quantum system 

2.1 Description of a particle in spin coordinates 

The basic equation of quantum mechanics is Schrodinger's equation, i.e. 

i^=H^(x,t) (I) 

where \ip(x,t)\ 2 is the probability density function of finding the particle at position x at time instant i, 
and H is the system's Hamiltonian, i.e. the sum of its kinetic and potential energy, which is given by 
H = p 2 /2m + V, with p being the momentum of the particle, m the mass and V an external potential. 
The solution of Eq. (JTJ) is given by ip(x, t) = e tHt ip(x, 0) II7I 

However, cartesian coordinates are not sufficient to describe the particle's behavior in a magnetic field 
and thus the spin variable taking values in SU(2) has been introduced. In that case the solution i\> of 
Schrodinger's equation can be represented in the basis \r, e > where r is the position vector and e is the 
spin's value which belongs in { — (fermion). Thus vector if) which appears in Schrodinger's equa- 
tion can be decomposed in the vector space |r, e > according to \tp >= Jd 3 r\r,e >,< r,e\i/j >. The 
projection of \ip > in the coordinates system r, e is denoted as < r, e\tp >— ijJe{r). Equivalently one has 
tp+{r) =< r, +\tp > and ip-(r) =< r,—\ip >. Thus one can write ip( r ) — [4 , + ( r )> - 0-( r )] T - 

2.2 Motion of the particle in a magnetic field 

The interaction between the particle's spin and a magnetic field and the effect this has on particles' motion 
is described in the Stern-Gerlach experiment. It can be observed that due to the gradient of the magnetic 
field and the particle's spin a deviation of the particle's trajectory takes places. It is assumed that the 
intensity of the magnetic field B z is positive while its gradient is negative. 

Since the particles are taken to be neutral they are not subject to Laplace forces F = q-vxB. The particles 
however have magnetic moment M which is associated to potential energy given by W — —M-B. There is 
also the kinetic moment F which is related to the magnetic moment M according to the relation M — j-T, 
where 7 is the gyromagnetic ratio. The force that is applied to the particle is the gradient of the particle's 
potential energy F = V(M-B). Between the kinetic moment T and the magnetic moment M the following 
relation also holds |£ = MxB=^^ = jTxB. The particle behaves like a gyroscope. The term ^ is 
perpendicular to F while the kinetic moment rotates round the magnetic field B z . It can be seen that the 
kinetic moment T is proportional and collinear to the magnetic moment M. According to F = \/{M-B) 
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the magnetic field B z forces the magnetic moment M to rotate round it, with constant angular velocity. 
To calculate force F the gradient of the particle's potential energy is found 

W = -M-B = -M X B X - M y B y - M Z B Z (2) 

It can be assumed that M z and M x = M y = 0. The latter ir true because the frequency of rotation of M 
is very high, and it is not possible for M x and M y to affect the particle's motion in any other way than 
their average value, which is 0. Then it holds 

F = V(M-B) = M Z VB Z (3) 

It also holds that = and = because the magnetic field is assumed to be independent of x and 
y. Therefore, the force which is responsible for the deviation HN (see Fig. Q]) of the particle from the 
straight line is proportional to the magnetic moment M z and is collinear to axis Oz. 

As the magnetic moments of the various particles are uniformly distributed between + |M| and — \M\ it 
is equally possible to find particles having magnetic moment between +\M\ and — \M\. Thus one expects 
that the particles beam will generate on plate P an equiprobable symmetric distribution diagram as shown 
with dashed line in Fig. [TJ However, in reality one observes two different distributions centered at points 
Ni and N2, which means that the particle's magnetic moment M z along axis z can take only two distinct 
values (magnetic moment eigenvalues or spin values). 




Figure 1: The two distributions denote that the magnetic moment of the particle can take only two distinct 
values (spin up or spin down) 



The results of the Stern-Gerlach experiment lead to the following conclusion: if some one measures the 
component of the kinetic moment T z of the particle (which is proportional to the magnetic moment M z ) 
then only one of the values corresponding to deviations HN\ and HN2 can be found. Therefore, the 
classical description of a magnetic moment vector M the measurements of which can take any value has 
to be abandoned. Consequently, the magnetic moment M z is a variable, with spectrum that can take only 
two eigenvalues (e.g. ±^). 
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2.3 Measurement operators in the spin state-space 



It has been shown that the eigenvalues of the particle's magnetic moment are ±i or iTij. The correspond- 
ing eigenvectors are denoted as |+ > and |— >. Then the relation between eigenvectors and eigenvalues 
is given by S z \ + >= +(h/2)\+ >, S z \— >= +(H/2)\— >, which means that the measurement of the 
Stern-Gerlach experiment shows only the two possible eigenvalues of the magnetic moment. In general the 
particle's state, with reference to the spin eigenvectors, is described by 

\tp >= a\+ > > 

(4) 

with \a\ 2 + \(3\ 2 = 1 
while matrix S z has the eigenvectors |+ >= [1, 0] and |— >= [0, 1] and is given by 



(5) 



_ h (i o 

Similarly, if one assumes components of magnetic moment along axes x and z, one obtains the other two 
measurement (Pauli) operators 



When the eigenvalue of the magnetic moment of the particle M z is at state |+ > or at state |— > the 
particle follows a well-defined trajectory as shown in Fig. [T] When the spin's state is a linear superposition 
of |+ > and |— > it is no longer possible to predict the trajectory of the particle. Therefore, there is a 
certain probability amplitude that the particle is found in one out of two trajectories. The probability to 
locate the particle at a specific point of surface P (particle's wave-function) has non-zero values at two 
different areas which are designated by points N\ and N 2 - Therefore, the particle may either appear round 
point Ni or round point iV 2 . However it is not possible to predict the particle's position with certainty, 
although the initial conditions are known (bifurcation). It is noted that in Fig. [T]there are not two different 
particles, but one particle, with a probability wave-function that consists of two parts centered at points 
Ni and N 2 . 



2.4 The spin eigenstates define a two-level quantum system 



The spin eigenstates correspond to two different energy levels. A neutral particle is considered in a mag- 
netic field of intensity B z . The particle's magnetic moment M and the associated kinetic moment T are 
collinear and are related to each-other through the relation M = ^T. The potential energy of the particle 
is W = —M Z B Z — —jB z T z . Variable luq = —-fB z is introduced, while parameter T z is substituted by the 
spin's measurement operator S z . 

Thus the Hamiltonian H which describes the evolution of the spin of the particle due to field B z becomes 
H = ojoS z , and the following relations between eigenvectors and eigenvalues are introduced: 

H\+>=+^\+>, JT|->=+£sa|-> (7) 
Therefore, one can distinguish 2 different energy levels (states of the quantum system) 

E + = +^|a, = -%f. (8) 

By applying an external magnetic field the probability of finding the particle's magnetic moment at one of 
the two eigenstates (spin up or down) can be changed. This can be observed for instance in the Nuclear 
Magnetic Resonance (NMR) model [17] , 



4 



3 Stochastic mechanics formulation of Schrodinger's equation 



3.1 Particles' motion can be written as a diffusion process 

Schrodinger's equation, which was described in Eq. (|TJ) , can be transformed into a diffusion equation by 
substituting variable it with t \V2\ . This change of variable results in the diffusion equation 

Eq. can be also written as ^ = —Hp, where H is the associated Hamiltonian and the solution is of 
the form p(x,t) = e~ p(x), and variable a 2 is a diffusion constant. The probability density function p 
satisfies also the Fokker-Planck partial differential equation |12| 



-± = {-y^~^)]p (io) 



dp 1 2 d 2 d 
~di~h a dx 2 ~dx~' 

where u(x) is the drift function, i.e. a function related through derivative to the external potential V. 
Now, as known from quantum mechanics, particle's probability density function p(x) is a wave-function 
for which holds p(x) — \ifj(x)\ 2 with tp( x ) = J2iLo c k' i l ; k(x), where ipk(x) are the associated eigenfunctions 
[T7] , It can be assumed that po(x) = \^ 2 {x)\, i.e. the p.d.f includes only the basic mode, while higher order 
modes are truncated, and the drift function u(x) of Eq. (TTU1) is taken to be u(x) — \o 2 df> g^ E2- 
Thus, it is considered that the initial probability density function is p(x) — po{x), which is independent 
of time. This means that the p.d.f. remains independent of time and the examined diffusion process is 
a stationary one, i.e. p(x,t) = po(x) Vt. A form of the probability density function for the stationary 
diffusion is that of shifted, partially overlapping Gaussians [12] . 



Continuing from Fokker-Planck's equation, given in Eq. (|10p . the Ornstein-Uhlenbeck diffusion is obtained 
which is a model of the Brownian motion [19] . The particle tries to return to the equilibrium x = under 
the influence of a linear force, i.e. there is a spring force applied to the particle as a result of the potential 
V{x). The corresponding phenomenon in quantum mechanics is that of the quantum harmonic oscillator 

(Q.H.O.) [IT] . Assuming a stationary p.d.f., i.e. p(x) — ipo(x) 2 = C 2 e~^, the force applied to the particle 
due to the harmonic potential (drift) V(x) is found to be u(x) = = —lux, which 

means that the drift is a spring force applied to the particle and which aims at leading it to an equilibrium 
position. 

Now, a kinematic model for the particles will be derived, in the form of Langevin's equation. The stochastic 
differential equation for the position of the particle is [12] : 

dx(t) = u(x(t))dt + dw(t) (11) 

where u(x) — —kx is the drift function, and is a spring force generated by the harmonic potential V{x) = 
kx 2 , which tries to bring the particle to the equilibrium x = 0. The term w(t) denotes a random force (due 
to interaction with other particles) and results in a Wiener walk. Knowing that the Q.H.O. model imposes 
to the particle the spring force u(x) = —ujx, Langevin's equation described in Eq. ()11[) . becomes 

dx(t) = -ujx(t)dt + dw(t) (12) 

Eq. (TT2l is a generalization of gradient algorithms based on the ordinary differential equation (O.D.E) 
concept, where the gradient algorithms are described as trajectories towards the equilibrium of an ordinary 
differential equation [18] . 
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3.2 Particle's spin in stochastic mechanics 



For a particle described in classical quantum mechanics by Schrodinger's equation, the spin represented 
the eigenvalues of the particle's magnetic moment (projection on the z axis) and the associated eigenstates 
were denoted as |+ > (spin up) and |— > (spin down). A representation for spin can be also obtained in 
terms of stochastic mechanics. The particle is described not only in the R 3 space by its cartesian coordi- 
nates, but is also described in the SU(2) space, i.e. it is described also by the spin variables which take 
two discrete values. 

As explained in the description of the Stern-Gerlach experiment in subsection 12. 2[ particles with spin ±2 
are emitted from a source and pass through a magnetic field B z . Then, the kinetic moment of the particle 
is proportional to its magnetic moment, and due to its interaction with the magnetic field generates a 
force which causes scattering of the particle's trajectory F — V(MB Z ). If the particle's path deviates 
towards the positive (negative) direction of the gradient then one can conclude that the particle's spin is 
at state |+ >, or |+ >. The particle's trajectory shows the discrete spin eigenvalues, while the particle's 
speed shows the discrete evergy levels E + and E- which correspond to the spin eigenvalues. It has been 
shown that starting from Schrodinger's equation and continuing to Focker-Planck and Ornstein-Uhlenbeck 
equations one obtains Langevin's stochastic differentinal equation, i.e. dxt = b(x,t)dt + dwt- It has been 
also shown that for long time t, b(x,t)~j, i.e. in b(x,t) — k-x one can consider k — i thus ones obtains 
the SDE [H] 

dxt — —x t dt + dw t (13) 

and based on the so-called "martingale convergence theorem" it has been proven that the limit lim t _- i . 00 
exists and is p + = lim^oo ^f. This limit is the kinetic moment which corresponds to the magnetic moment 
eigenstate |+ > and to the magnetic moment eigenvalue 'spin- up'. Moreover, it has been shown that for 
every measurable subset Aei? 3 , the probability P+(pt £ A) is equal to the quantum mechanical probability 
that the final moment of the particle belongs in A. Thus, in stochastic mechanics a way to measure the 
moment is through the limit p + = lim^oo^. 

Consequently, in stochastic mechanics, the different energy eigenstates of the particle can be conceived 
according to the Stern-Gerlach experiment as follows: the particle has initial moment p_, and while it 
approaches to the surface P in which points Ni and JV2 belong (see Fig. [IJ, the particle's trajectory 
becomes straight and the final moment becomes p+. Then, the only possible values for the change of 
energy are 

m(| P+ | 2 -b_| 2 ), i.e E+-E- (14) 

where, Ej is the Hamiltonian's eigenvalue, as explained in Eq. [5] Therefore, in stochastic mechanics the 
stage of the particle's motion at which its trajectory becomes straight and its moment stabilizes at the 
final value p± can be considered as a collapse of the particle's wave function. 

4 Open-loop control scheme for a multi-particle system 

As explained in subsection 12 . 2 1 the change in the particle's moment means transition between the particle's 

energy eigenstates. This can be a manner to control the probability to find the particle's magnetic moment 

in one of the spin's eigenstates. The particle's kinematic model has been described by Eq. (fl2"|) . and here 

it will be formulated as follows: the motion of the particle is described by the stochastic oscillator model 
,2 1 

TO ^ti + c^" = /, where x is the position coordinate, m is the mass, c is the coefficient of viscous friction, 
and / is the aggregate force acting on the particle [2Uj . Defining x\ = x and X2 — particle's motion can 
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be written in state-space form = x 2 , and ^ = — + hf- Keeping the second of the state-space 
equations one has the Langevin equation for the particle's velocity 

% = --v+ -f. (15) 

at m m J v ' 

Next, an open-loop control, based on flatness-based control theory, will be applied to the particles: 

Definition: The system x — f(x,u), x <E R n ,u € R m is differentially flat if there exist relations h : 
R n x(R m ) r+1 ^R m , (j) : {R m ) r ^R n and ip : {R m ) r+1 -» R m , such that y = h(x,u,u,--- x = 

(f>(y, y, ■ ■ ■ , y^ -1 ') and u = ip(y, V, ' ' ' 1 2/ ; 2/ )■ This means that all system dynamics can be expressed 
as a function of the flat output and its derivatives, therefore the state vector x and the control input u can 
be written as x(t) = <P(y(t),y(t), • • • , yW(t)) and u(t) = i>(y(t),y(t), • ■ ■ , y^ r+1 \t)) UUMM- 
If a control term u 1 is introduced in Eq. (|15[) this can be written as: 

v l = —lov 1 + u l + if (f 6) 

where —ujx % is the drift term due to an external potential, u l is the external control and rf is a disturbance 
term due to interaction with the rest iV-1 particles, or due to the existence of noise. Then it can be easily 
shown that the system of Eq. (| 16|) is differentially flat, while an appropriate flat output can be chosen to be 
y = v l . Indeed all system variables, i.e. the elements of the state vector and the control input can be writ- 
ten as functions of the flat output y 7 and thus the model that describes the z-th particle is differentially flat. 

A control input that makes the i-th particle track the reference trajectory y\ is given by 

u L = ujv l r + «* + u*, (17) 

where v l r is the reference velocity profile for the i-th particle, and v\ is the derivative of the i-th desirable ve- 
locity. Moreover u l c = —if stands for an additional control term which compensates for the effect of the noise 
r\i on the i-th particle. Thus, if the disturbance r\i that affects the ith-particle is adequately approximated 
it suffices to set u % c = —rji. The application of the control law of Eq. (jTTJ) to the model of Eq. (|16|) results in 
the error dynamics v l — v], — uiv' 1 + ujv], + rf — it* =h i- e - v l — v\ +ui(vi — vl) = rf + u c =>■ e l +uie l = rji + u c . 
Thus, if u c = —rji then limt^ooe = 0. 

Next, the case of the N interacting particles will be examined. The control law that makes the mean of 
the multi-particle system follow a desirable velocity profile E{v l r } can be derived. The kinematic model of 
the mean of the multi-particle system is given by 

E{v 1 } = -uE{v'} + E{u 1 } + Elf} (18) 

i = 1, • ■ • ,N, where E{v 1 } is the mean value of the particles velocity, E{v 1 } is the mean acceleration, E{rf} 
is the average of the disturbance signal and E{u 1 } is the mean control input. The open-loop controller is 
selected as: 

£{y} = uE{v% + E{v% - E{rf) (19) 

where E{v l } r is the desirable mean velocity. Assuming that for the mean of the particles system holds 
E{rf} = 0, then the control law of Eq. (|19p results in the error dynamics E{e 1 } + uiE{e 1 } — 0, which 
assures that the mean particles' velocity will track the desirable velocity profile, i.e. \\m t ^, oE{e' 1 } = 0. 
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5 Conclusions 



The paper has proposed the spin-orbit interaction as a model of quantum computation. The stochastic 
mechanics formulation of Schrodinger's equation was introduced. This enables to consider the control of 
the entangled spin state as a problem of motion control of the mean of a particles ensemble along a reference 
trajectory. It is shown that using an appropriate open-loop control (which can be created by a magnetic 
field) the mean velocity of the particles can track a desirable velocity profile. This in turn implies that the 
quantum state represented by the mean particles moment can be also controlled. 
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